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We study the angular structure of polarization of light transmitted through a ne- 
matic liquid crystal (NLC) cell by analyzing the polarization state as a function of 
the incidence angles and the polarization of the incident wave. The polarization re- 
solved angular (conoscopic) patterns emerging after the NLC cell illuminated by the 
convergent light beam are described in terms of the polarization singularities such 
as C-points (points of circular polarization) and L-lines (lines of linear polarization) . 
For the homeotropically aligned cell, the Stokes polarimetry technique is used to 
measure the polarization resolved conoscopic patterns at different values of the ellip- 
ticity of the incident light, eg™^^ impinging onto the cell. Using the exact analytical 
expressions for the transfer matrix we show that variations of the ellipticity, Cg^'^'', 
induce transformations of the angular pattern exhibiting the effect of avoided L-line 
crossings and characterized by topological events such as creation and annihilation 
of the C-points. The predictions of the theory are found to be in good agreement 
with the experimental results. 
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I. INTRODUCTION 



Singularities that represent structurally stable topological defects have long been known 
to play a pivotal role in condensed matter physics [T]. They are of particular importance in 
determining the properties of systems with a broken continuous symmetry [21 El IH E] (for a 
recent review see, e.g., Ref. [B]). Topological methods of quantum field and gauge theories 
have been extensively used to classify the topological defects and to describe transformations 
of the singularities in ordered media such as superfiuids and liquid crystals [Zl El |9]. 

For electromagnetic vector fields, it was originally recognized by Nye [101 [HI [12] that 
the so-called polarization singularities are the important elements characterizing geometry 
of the Stokes parameter fields. In particular, the polarization singularities such as the C- 
points [the points where the light wave is circularly polarized] and the L-lines [the curves 
along which the polarization is linear] frequently emerge as the characteristic feature of 
certain polarization state distributions. Over the past two decades these singularities and 
related issues have been the subject of numerous theoretical, experimental and numerical 
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studies [l3l[Tl[l5l[T6l[I71[T8l[T9l|2ni[2ll|22l[^ 

The theory of polarization singularities has also been found to be a useful tool for study- 
ing optical properties of anisotropic media. In Ref. ^S], it was applied to study the an- 
gular dependence of the polarization state of the electric displacement field for plane wave 
eigenmodes in birefringent dichroic chiral crystals. This analysis was then generalized and 
extended to a more complicated case of bianisotropic media |26] . 

The experimental results and theoretical analysis presented in Refs. [271 EH] deal with 
the unfolding of a linearly polarized Laguerre-Gauss (LGqi) beam with an on-axis vortex 
on propagation through a birefringent crystal. It was found that a complicated pattern of 
polarization singularities is formed as a result of the anisotropy induced symmetry breaking. 

In our recent investigation into the two-dimensional angular distributions of the Stokes 
parameters describing the polarization structure behind the conoscopic images we performed 
a comprehensive theoretical analysis of the polarization state of the light transmitted through 
nematic liquid crystal (NLC) cells as a function of the incidence angles [291 130] . This 
polarization structure — the so-called polarization resolved angular ( conoscopic ) pattern — 
is represented by the field of polarization ellipses and results from the interference of four 
eigenmodes excited in NLC cells by the plane waves with varying direction of incidence. 

Note that NLCs are technologically important materials where the optical anisotropy 
is determined by the orientational structure which is sensitive to external fields and, in 
restricted geometries, can also be influenced by changing the boundary conditions [311 E21 
[55] . The conoscopy is widely used as an experimental technique to characterize orientational 
structures in NLC cells. 

For example, this method was employed to detect biaxiality of NLCs [Ml EH] and to 
measure the pretilt angle in uniaxial liquid crystal cells [361 EZ]- Orientational structures 
and helix unwinding process in ferroelectric and antiferroelectric smectic liquid crystals were 
also studied by conoscopy in Refs. [381 EHl HQ]. 

In this paper we explore both theoretically and experimentally the polarization resolved 
angular patterns for convergent light beam impinging on the homeotropically aligned NLC 
cell at varying polarization of the incident wave. The layout of the paper is as follows. 

In Sec. [IT] we give experimental details and describe our setup employed to carry out 
measurements using a suitably modified method of the Stokes polarimetry [271 SIl 1121 113] ■ 

The problem of light transmission through a uniformly anisotropic NLC cell is considered 
in Sec. Ill A Using the 4x4 matrix formalism we deduce the general expressions relating 
the evolution matrix to the transmission and reflection matrices. The analytical results are 



used for analysis of the polarization resolved conoscopic patterns, in Sec. Ill B 



The patterns emerging after homeotropically aligned NLC cells are treated in Sec. Ill C 
The polarization ellipse fields are studied as a function of the ellipticity, el^^^\ and the 

polarization azimuth, 4'p^'^\ characterizing the polarization state of the incident wave. 

The analytical expressions describing loci of the C-points and L-lines along with the 
index formula are used to examine rearrangements of the polarization singularities caused 
by variations of the polarization parameters e^J^'^^ and (f)p^'^\ 

We find that changing the azimuth results in rotation of the ellipse field as a whole by the 
angle (f)p^^\ whereas the transformations induced by the ellipticity appear to be complicated 
by the presence of bifurcations leading to creation and annihilation of the C-points. It is 
shown that the structure of the intersecting L-lines formed at e^^ '^^ = smoothly evolve into 

a family of the concentric L-circles at |eeii"^'*| = 1. At small values of the ellipticity, the effect 
of avoided L-line crossings is found to be an important feature of this transition. 
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Experimentally measured and computed fields of the polarization ellipses are presented in 
Sec. |IV| along with discussion and concluding remarks. Mathematical details on the evolution 
operator of uniformly anisotropic media and related issues are relegated to Appendix |Aj 

L,PH Li P.Wj Lg WjJS L^^^^a^ 

FIG. 1: Experimental setup: He-Ne is the laser; Li, L2 and L5 are the collimating lenses; L3 
and L4 are the microscope objectives; PH is the pinhole; NLC is the NLC cell; Pi and P2 are 
the polarizers; Wi and W2 are the quarter wave plates; CCD is the CCD camera. A microscope 
objective L3 is illuminated with an elliptically polarized and expanded parallel beam of light from 
a He-Ne laser. The output from a second objective L4 is collected by a CCD camera through the 
Stokes analyzer. 



II. EXPERIMENT 



In our experiments we used the NLC cells of thickness d = 110 /xmfilled with the nematic 
liquid crystal mixture E7 from Merk. Two glass substrates were assembled to form a home- 
tropically oriented NLC cell. At the wavelength of light generated by a low power He-Ne 
laser from Coherent Group with A = 632.8 nm (see Fig. [T]), the ordinary and extraordinary 
refractive indices of the NLC are Uo = 1.5246 and Ue = 1.7608, respectively; the refractive 
index of the glass substrates is Ug = 1.5. 

Figure [T] shows our experimental setup devised to perform the conoscopic measurements 
using the Stokes polarimetry technique SU 1121 SS] • Referring to Fig. [T| the cell is irra- 
diated with a convergent light beam formed by the microscope objective of high numerical 
aperture L3. The input polarizer Pi is combined with the properly oriented quarter wave 
plate Wi to control the polarization characteristics (the ellipticity and the azimuth of polar- 
ization) of a He-Ne laser beam which is expanded and coUimated using the lenses Li and L2. 
A charge coupled device (CCD) camera collects the output from the microscope objective 
L4 through the collimating lense L5 and the Stokes analyzer represented by the combination 
of the quarter wave plate W2 and the polarizer P2. 

This optical arrangement therefore collects simultaneously the transmittance of the cell 
for a range of incident angles. The distribution of the Stokes parameters 5*0 ... 5*3 describing 
the state of polarization of the transmitted light can then be obtained by performing the 
measurements at six different combinations of the quarter wave plate W2 and the analyzer 
P2 and using the well-known relations [44] 

5*1 = Iqo — /900, (la) 

S2 = -^45° " h35°y (lb) 
>S'3 = /rcP — -^LCP, (Ic) 

So = ^Sl + SI + SI (Id) 
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where Iqo , /ggo , /450 and /1350 are the hnearly polarized components (subscript indicates the 
orientation angle of the analyzer); /rcp and /lcp are the right- and left-handed circular 
polarized components measured in the presence of the quarter wave plate W2. 

Geometrically, the distribution of the Stokes parameters measured in the observation 
plane can be conveniently represented by the two-dimensional field of the polarization el- 
lipses. The geometrical elements of polarization ellipses that determine the polarization 
ellipse field can be readily computed from the Stokes parameters ([T]). 

The orientation of a polarization ellipse is specified by the azimuthal angle of polarization 
{polarization azimuth) 

0p = ^ arg(5i + 182) (2) 
and its eccentricity is described by the signed ellipticity parameter 



eeii = tan 



1 . fS3 



(3) 



This parameter will be referred to as the ellipticity. The handedness of the ellipse is deter- 
mined by the sign of the ellipticity parameter eeu. 



In Sec. IV, we shall describe the polarization ellipse fields representing the experimentally 
measured angular patterns. These will be compared with the predictions of the theory 
discussed in the subsequent section. 



III. THEORY 

A. Transmission boundary-value problem 

We consider a nematic liquid crystal (NLC) cell of thickness d sandwiched between two 
parallel plates that are normal to the z axis: z = and z = d. Typically, anisotropy of nemat- 
ics is locally uniaxial and NLC molecules align on average along a local unit director, d [31] . 
In this case the NLC director 

d = eo(d) = sin 6'd cos 0d x + sin 9^ sin 0d y + cos 6'd z (4) 

determines the optic axis. So, the uniaxially anisotropic dielectric tensor, e, is given by 

s = e_Ll3 + Aed O d, Ae = - e_L, (5) 

where I„ is the n x n identity matrix. (In what follows carets will denote unit vectors.) 
Its two principal values e± and ey define the ordinary and extraordinary refractive indices, 

= -^/ie± and Ue = ^/JJ^\, where /i is the NLC magnetic permeability. 

In a more general case of biaxial nematics [13 HB] that were recently observed experi- 
mentally [SnillTj, there are three different dielectric constants ei, €2 and €3 representing the 
eigenvalues of the dielectric tensor: 

£ = €3 13 + Aei d (g) d + Ae2 m ® m, (6) 

where Ae^ = — €3, and the eigenvectors d, m and 1 = d x m give the corresponding 
principal axes. 
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Similar to the director Q, the unit vectors m and 1 can be expressed in terms of Euler 
angles as follows 

m = cos 7d ei (d) + sin 7d 62 (d) , (7a) 
i = -sin7dei(d) + cos7de2(d), (7b) 

where ei(d) = (cos 6'd cosf^d, cos^d sin(/)d, — sin 6'd) and 62(0!) = (— sin(/)d, cos(/)d, 0). 




FIG. 2: Geometry of nematic cell in the plane of incidence. 



We shall need to write the Maxwell equations for a harmonic electromagnetic wave (time- 
dependent factor is exp{—iujt}) in the form: 



V X E = i/ifcvacH 

V X H = 



(8a) 
(8b) 



where /Cvac = uj/c is the free-space wave number; /j, is the magnetic permittivity and D = e-E 
is the electric displacement field. 

The medium surrounding the NLC cell is assumed to be optically isotropic and char- 
acterized by the dielectric constant em and the magnetic permittivity /im. So, Maxwell's 
equations in the medium outside the cell can be obtained from Eq. ([s]) by replacing /i and 
D with and emE, respectively. 
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As is shown in Fig. |2| there are two plane waves in the half space z < bounded by the 
entrance face of the NLC cell: the incoming incident wave {Ejnc, Hinc} and the outgoing 
reflected wave {Erefi, Hrefl}. In the half space z > d after the exit face, the only wave is the 
transmitted plane wave {Etr,Htr} which propagates along the direction of incidence and is 
excited by the incident light. 

So, the electric field outside the cell is a superposition of the plane wave solutions of the 
Maxwell equations 

EUo = Ei„e(kinc)e'('^''^-^) + E,efl(krefi)e^(^-«-'-), (9a) 
EU, = Etr(ktr)e*(''-^\ (9b) 

where the wave vectors kinc, krefi and ktr that are constrained to lie in the plane of incidence 
due to the boundary conditions requiring the tangential components of the electric and 
magnetic fields to be continuous at the boundary surfaces. These conditions are given by 

P(z) ■ [E|,=o+o - E|,=o-o] = P(z) ■ [E|,=,+o - E|,=,_o] = 0, (10a) 
P(z) ■ [H|,=o+o - H|,=o-o] = P(z) ■ [H|,=,+o - H|,=rf_o] = 0, (10b) 

where P(z) = I3 — z z is the projector onto the plane with the normal directed along the 
vector z (the x-y plane). 

Another consequence of the boundary conditions ( 10 ) is that the tangential components 
of the wave vectors are the same. Assuming that the incidence plane is the x-z plane we 
have 

ka = /CvacQa = ^mka = fc^ X + k'f^ z, « G {iuc, refl, tr}, (11) 

where fem/Z^vac = = i/A^mem is the refractive index of the ambient medium and the 
components can be expressed in terms of the incidence angle 6*^0 as follows 

kx = sill ^inc = k^ncQxi (12) 
= ^2 = k^ = fcjn COS ^inc = k^gjc Q^z ^- (l^) 



The plane wave traveling in the isotropic ambient medium along the wave vector (11) is 
transverse, so that the polarization vector is given by 

E„(kJ = Ej")ei(k„) + E^^^e^iK), (14) 
ei(k,) = A;^i(A;(-)x-A;,z), e2(kj = y, (15) 

where E^^'^ and -E^"'' are the in-plane and out-of-plane components of the electric field, 
respectively. The vector characterizing the magnetic field is 

/x^ H„(k«) = X E,(k,) = [e\^^ - El")ei(k„)] , (16) 

where = k~^^a = 'T-mka- Note that, for plane waves, the dimensionless vector 

q = fc^ack (17) 

is parallel to k and its length gives the refractive index. For convenience, we shall often use 
this vector in place of the wave vector. 
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The expressions ([9])-(16) give the electromagnetic field of incident, transmitted and re- 
fiected waves propagating in the ambient medium. This field is of the general form 



{E(r),H(r)} = {E(z), H(^)} exp[2(kp, r)], kp 

qP = ^vac kp = qp [cos(0inc) X + sin(0inc) f] , 



Pfz)-k, 



(18) 
(19) 



where the azimuthal angle 0inc specifies orientation of the incidence plane. So, in our case, 
we have 0inc = and qp = q^- 

The representation (18) can be applied to describe the field inside the cell when the 
dielectric tensor is independent of the in-plane coordinates x and y. In this case the lateral 
components of the electric and magnetic fields 



Ep = P(z)-E = E,x + E,y = , 



Hp = H X z = Hy± - H^y 
can be conveniently combined into the vector 



Hy 



(20) 
(21) 



Fiz) 



Ep{z) 
Hp(^) 



Ey{z) 

Hy{z) 

\-Hx{z)) 



(22) 



From Maxwell's equations ([8]), the components normal to the substrates, and Hz, can be 
readily expressed in terms of the lateral components 



^iHz = (qp, Ep X z) , e^zEz = - (z, £ ■ Ep) - (qp. Hp) . 



(23) 



By using the relations (|23|) to eliminate the normal components we obtain the equation for 

(24) 



the field vector (22) in the following matrix form 
-idrF = M • F E 



Mn Mi2 
M21 M22 



E p 
Hp 



T — ^vac^^! 



where 



M. 



(22) 



-zz ^ctz 



(12) 



af3 



'^zz qa qp 1 



M. 



(21) 



a/3 



-1 (P) (P) -1 



(25a) 
(25b) 
(25c) 
(25d) 



where Sap is the Kroneker symbol, qp = k^^^'kp and pp = z x qp. 



The field vector (22) can now be written in the form of a general solution to the linear 

F(r) = U(r) ■ F(0), (26) 



problem (24) 
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where U(r) is the evolution operator that can be determined by solving the initial value 
problem 

-idrV = M ■ U, U(0) = I4. (27) 

Our above notations can be regarded as a version of the well known matrix formalism 
due to Berreman jH]. By using these notations we derive the boundary conditions (10) in 
the matrix form 



F< 
F> 



E- 

-'-'mc 

Ercfl 





= F(0), 

F{h), h = Ks.cd, 



(28a) 
(28b) 



where E„ = 



E 



The matrix Vm relates the field vectors, F< = F(0 — 0) and F> = F{h + 0), and the 
vector amplitudes E^ of the waves in the surrounding medium. From Eqs. (14)- (16), its 
block structure 



E„ 



— CTsEm 



is characterized by the two diagonal 2x2 matrices 

Em = diag(g("'V'^m, 1), Aim Hm = diag(nm, q'f'^] 
where gi™"* = ^/n"^ — — diag(l, —1). 



(29) 



(30) 



trix (29) 



It is rather starightforward to check the validity of the algebraic identity for the ma- 

(31) 



G V, 



A/'mdiag(l2, -I2) 



G = 



I2 
I2 



where N,^ = 2gi™''//Xni and the superscript T indicates matrix transposition. 

According to Ref. |19], identities of the form (31) can be derived as the orthogonality 
relations resulted from the conservation law for the energy flux in non-absorbing media. 
Algebraically, Eq. (31) can be used to simplify inversion of the matrix 

Vr„' = ■ [V„/ • G (32) 

and to ease qualitative analysis [see , e.g., Ref. for a more extended discussion of appli- 
cations] . 



After substituting Eq. (|26|) into the boundary conditions (|28a|) we have 

= W 



E- 

-'-'mc 

, Erefl 

w 



Etr 





(33) 
(34) 
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where U'^ = V-\h). 
From Eq. (33) we have 



E 



(try 
ftr) 




(35) 



(36) 

where T is the transmission (transfer) matrix hnking the transmitted and incident waves. 
Similar result for the reflected wave is 



^1 
E 



(refl)' 
refi) 



R 



E 



(inc)' 
(inc) 



R = W21 • Wr/ = W21 ■ T, 



(37) 
(38) 



where R is the reflection matrix. 



For uniformly anisotropic media, the elements of the matrix (25 ) are constants and solving 



Eq. (27) gives the formula for the evolution matrix 

U(r) = exp{iMr} = V ■ exp{iAr} ■ V~\ 
A = diag(Ai, A2, A3, A4), M ■ V = V ■ A, 



(39) 
(40) 



where Aj is the eigenvalue of the matrix M and V = (Vi V2 V3 V4) is the matrix composed 
of the right eigenvectors, M ■ Vj = Aj Vj. Technical details on computing the evolution 
operator for both biaxial and uniaxial anisotropies can be found in Appendix |X} In this 
section we restrict ourselves to the simplest (and important for our purposes) case that 
occurs when the NLC cell is homeotropically aligned. 

For the homeotropic director structure with d = z, the eigenvalue problem for the matrix 
M is easy to solve. The result is 

El — cTsE I 



A 



diag(g(^),gi°),-g(^),-g. 



diag(g^''V^^±, 1), AtH+ = diag(?2_L, qf^) 



V 



H 



(oh 



(41) 
(42) 



where qi"^ = ^/n^ — q^ and q. 



2 



- qi- 



The analytical relations (41) can now be substituted into the expression (39) to yield the 
evolution matrix \J~^{h) = \J{—h) that enter the linking matrix W defined in Eq. (34). 
The final step involves substituting the block 2x2 matrices Wn and W21 into Eqs. (36) 
and (38) to find the transfer and the transmission matrices in the following form 

diag(te,to) 



T 
R 



where 



diag(re, ro] 
A+/i = diag(g(^),g(°))/i = diag(5„5„ 

r± = (r±r-i)/2. 



cos(A+/i) — zr_|_ ■ sin(A_|_/;,) 
-zcr3 ■ r_ ■ sin(A+/i) ■ T, 



/^2„(m) 



2 (e) 



V 









q^\ 



= diag(7e,7o). 



(43) 
(44) 

(45) 
(46) 
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3 (b) Lemon: Ic = +1/2 and 
Nc = 1 



(a) Star: Ic = -1/2 and Nc 
FIG. 3: Arrangment of the polarization ellipses around the C-points of three different types. 



(c) Monstar: Ic = +1/2 and 
Nc = 3 



B. Polarization resolved angular patterns 

We will now proceed to a study of how the polarization properties of the transmitted wave 
depend on the direction of the incident beam. This direction is specified by two angles: the 
incidence angle 6'inc and the azimuthal angle of the plane of incidence 0inc- 

In the previous section and in Appendix |A| the transmission problem was analyzed in 
the plane of incidence where 0inc = 0. Clearly, when (p^^^ 7^ 0, in order to have the Euler 
angles describing orientation of the director Q with respect to the incidence plane we need 
to replace the director azimuthal angle 0d with 0d — (p^^c- 

In this section dependence of the polarization parameters of transmitted waves on the 
angles 6'inc and 0inc will be of our primary concern. For this purpose, we assume that the 



transmission matrix T(6'inc, Od, 0d) considered in Sec. Ill A and in Appendix [A| is changed to 



T(6'inc, 6'd, 



As previously discussed in Sec. [IT], experimentally, the characteristics of the polarization 
ellipse can be determined by measuring the Stokes parameters ([T]). These parameters are 
related to the coherence matrix with the elements Mq,/3 = EaEj^*, where a,f3 G {||,-L 
} [5T| [52] . In the circular basis, \/2e±(k) = ei(k) ± ie2(k), this matrix written as a linear 
combination of the Pauli matrices gives the Stokes parameters. Si, as its coefficients 



C M C+ 



\E+\^ E+El 
E^El \E\^ 



(47) 



where 



(To = I 



2, 



E_ 



C 



1 

1 



CT2 



E\\ 
E, 



C 



^1/2 



-i 

1 

1 - 

1 ? 



0-3 



1 
-1 ' ' 



(48) 
(49) 



Since the determinant of the coherence matrix vanishes, det M = 0, the Stokes parameters 
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lie on the four-dimensional cone Sq = Yl^=i ^'^d can be parameterized as follows 

Sq = \E+\^ + |E„|2 = + \E^\^, (50a) 

51 = 2ReElE_ = j^nl^ - = 5ocos2xpCos20p, (50b) 

52 = 2 Im = 2 Re E^E^^ = Sq cos 2xp sin 20p, (50c) 
^3 = - = 2ImE^E|| = S'osin2xp, (50d) 

where < cpp < ir is the polarization azimuth ^ and — 7r/4 < Xp < 7r/4 is the elliptic- 
ity angle. Then, the relations expressing the ellipse characteristics in terms of the Stokes 
parameters are 

0p = 2-1 arg(E;^_) = 2-1 arg^, S = Si + 1S2, (51) 
eeii = , I ~| = tanxp, Xp = 2-^ arcsin(S'3/S'o). (52) 



Similar to Refs. [211 ES], we have used Eq. (47) and Eq. ( |50d[ ) to define the Stokes 
parameter S3 which is opposite in sign to that given in the book [3^ . The ellipse is considered 
to be right-handed (RH) if its helicity is positive, so that \E+\ > \E^\ and eeii > 0. In the 
opposite case with eeii < 0, the ellipse is left-handed (LH). 



In the spherical basis, {e+,e_}, similar to Eq. (35), the transmission matrix 



C ■ T • C+ (53) 



relates the circular components of the incident and transmitted waves 



E^^'n~ ^ \e 



M I = Tc ■ J„c) 1 , (54) 



where E^^ = (-E'||"'* =F iE^f^) / \/2. So, for the transfer matrix describing the conoscopic 
patterns on the transverse plane of projection, we have [291 EQ] 



T(p,0) = exp(-i0cr3) ■ Tc(p,0) ■ exp(i0cr3), (55) 

p = r tan 6'inc, = 0inc, (56) 

where p and are the polar coordinates in the observation plane {x = p cos and y = p sin 
are the Cartesian coordinates) and r is the aperture dependent scale factor. 

For the elliptically polarized incident plane wave with the circular components 

4'-) = exp(-^z.0f ^)) [1 + z/eiir^] li^incl, (57) 

expressed in terms of the polarization azimuth 0p"'^^ and the ellipticity parameter e^^^^\ the 
reduced components of the transmitted wave are given by 

t,,_4l-z/ei;r^]exp(-22z/V')}exp(-^z/0f=)), (58) 
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,{tr) 



where ip = (f) — (j)p^^^ . 

For the transmitted wave, the formula (58) gives the normahzed Stokes parameters s\ 
S^^^^ / Sq'^^ as a function of the incidence angles, 6i^c and 0inc- From Eqs. (51 ) and (52), these 
parameters are determined by the ellipticity, eeu = — + |^-|), and the 

polarization azimuth, (pp = — 0+)/2. Owing to the relation (56), the incidence angles and 



the points in the observation plane are in one-to-one correspondence. So, the distribution 
of the normalized Stokes parameters can be evaluated by computing the characteristics of 
the polarization ellipse eeii and (j)p at each point of the projection plane. Geometrically, this 
procedure yields the two-dimensional field of polarization ellipses that might be called the 
polarization resolved angular (conoscopic) pattern. 

The point where = and thus the transmitted wave is circularly polarized with 
eeii = — z/ will be referred to as the C^-point. This is an example of the polarization singularity 
where the phases (pu and (j)p become indeterminate. 

Equivalently, C-points can be viewed as the phase singularities of the complex scalar field 



S = 



Si + 182 



(59) 



proportional to the Stokes field defined in Eq. (51). Such singularities are characterized by 



the winding number which is the signed number of rotations of the two-component field 
{Si, 82) around the circuit surrounding the singularity [2J. The winding number also known 
as the signed strength of the dislocation is generically ±1. 

Since the polarization azimuth (51) is defined modulo vr and 20p = argS", the dislocation 
strength is twice the index of the corresponding C,y-point, Ic- For generic C points, Ic = 
±1/2 and the index can be computed from the formula 



sign 



lm{dJ*dyS) 



X=Xi, ' 



(60) 



where d^f is the partial derivative of / with respect to x. 



The relation (60) gives the index of the Cy-point with the coordinates {xuiUu) expressed 
in terms of the vorticity [iHl El]: hii{dxS*dyS) = dxSidyS2 — dySidxS2. The formula 
linking gradients of the complex field (E, E) oc \1/4.\1/_ and the index for C-lines in the 
three-dimensional space was derived in Ref. [54j. 



In the case with = (Ci,-point), only derivatives of '^y enter the expression (60) which 



can be suitably rearranged to yield the index of the C,y-point in the following form: 



sign 



\m{dx^,dy^!l) 



V 



sign 



y=yv 



p=Pu 



(61) 



Subsequently, we shall apply the formula (61 ) expressing the index in terms of the derivatives 



with respect to polar coordinates to the case of homeotropically aligned cell. 

In addition to the handedness and the index, the C-points are classified according to the 
number of straight lines terminating on the singularity. This is the so-called line classification 
that was initially studied in the context of umbilic points 

For generic C-points, as is shown in Fig. [3} the number of the straight lines, Nc, may 
either be 1 or 3. This number is 3 provided the index equals —1/2, Ic = —1/2, and such 
C-points are called stars. At Ic = 1/2, there are two characteristic patterns of polarization 
ellipses around a C point: (a) lemon with Nc = 1 and (b) monstar with Nc = 3 [10]. 
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FIG. 4: The real part of the function (69) that enter the left-hand side of Eq. (68) computed for 
(a) the cell filled with the NLC mixture E7 and (b) the KDP (potassium dihydrogen phosphate) 
birefringent crystal. The parameters are: (a) n± = 1.5246; n|| = 1.7608; = 1.5; d = 110 fim; 
and (b) n± = 1.5093; nn = 1.4683; = 1-0; d = 300 /im. 



Different quantitative criteria to distinguish between the C-points of the lemon and the 
monstar types were deduced in Refs. [211 15U] . 

The location of the C,y-points on the projection plane is determined by the polar coor- 
dinates: p^i^^ and 4'^^\ where k is the numbering label. These can be found by solving the 
equation 

|vl>,(p,0)| = O (62) 

that generally has multiple solutions. 

The case of linearly polarized wave with egu = provides another example of the polariza- 
tion singularity where the handedness is undefined. The curves along which the polarization 
is linear are called the L-lines. 

The transmitted wave is linearly polarized when the condition 

|M/+(p,0)| = |M/_(p,0)| (63) 



is satisfied. So, Eq. (63) describes loci of points forming the L-lines lying in the projection 
plane. 
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(b) 



FIG. 5: Schematic respresentation of creation and annihilation of C-points governed by saddle- 



node bifurcations of the solutions of Eq. (68) that take place within the areas marked by the 
boxes in Figs.|4j^a) and|4]j^b). (a) C-points annihilate as the governing parameter /3 passes through 
its critical value /3c in the immediate vicinity of unity, (b) Two successive bifurcations result in 
creation and annihilation of C-points at (3 = Pb and P = (3u, respectively. 



C. EUipticity induced effects 



1. C-points: creation and annihilation 
When the director is normal to the substrates, the NLC cell is homeotropically aligned 



and dz = 1. In this case the transmission matrix (43) is diagonal and its elements in the 
circular basis are given by 



t^, t± = (te±to)/2, (64) 

where the transmission coefficients 

te,o = te,o{qx) = [cOs{6e,o) - -116,0 Sm{6e,o)] ^=te,o{p) (65) 

are expressed in terms of the phases 



5e = q. 



-qlh 



(66) 



and the amplitudes 7e,o defined in Eq. (46). 



Owing to the cylindrical symmetry of the homeotropic structure, the transmission coef- 

E 0inc. From 



ficients (65) do not depend on the azimuthal angle of the incidence plane 



Eq. (58) another consequence of this symmetry is that the azimuth. 



t{inc) 



and the 

ellipticity of the transmitted wave depend only on the difference of the azimuthal angles: 
. It follows that the sole effect of changing the polarization azimuth of the incident 
wave is the rotation of the polarization ellipse field by the angle 



/,(inc) 



So, we focus our 



attention on the effects governed by the ellipticity of the incident wave, 



From Eqs. (58) and (64), the Ci,-points may appear only if the condition 



\t\ 



1 I (inc) 

1 + K\\ 

(inc) 



(67) 



cll 
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FIG. 6: Polarization resolved conoscopic patterns computed as polarisation ellipse fields in the 
observation plane for the homeotropically oriented cell filled with the NLC mixture E7 at small 
values of the ellipticity: (a) Cg™"^^ = 0.0, (b) e^™"^^ = 0.01, and (c) e^^^^^ = —0.01. The parameters 
used in calculations are: the wavelength is Ainc = 632.8 nm; the polarization azimuth is (j)p^^^ = 0; 
the cell thickness is d = 110 ^m; the scale factor is r = 5 mm; = 1.5, n± = 1.5246, and 
ny = 1.7608. The C-points of the star, the lemon and the monstar types are marked by circles, 
diamonds and squares, respectively. L-lines are represented by solid lines. Left-handed and right- 
handed polarization is respectively indicated by solid and open ellipses. 



is satisfied. Another form of this condition is 



2e 



(inc) 
cU 



Re F{p) = -V f--- = eeii /ii„c /?, (68) 

F(n) = ^'^^P)'*o^P) = ^^P(^^) (69) 

\Up)\^ + UpW cosh(ln|te|-ln|to|)' ^ ^ 
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FIG. 7: Polarization resolved conoscopic patterns of the homeotropically oriented NLC cell com- 



puted at different values of the ellipticity. Three cases are shown: (a) e. 



< e 



< e 



Fig.|6} 



(inc) 
ell 



< e 



1 ' 



(c) 



< e 



(inc) 
ell 



(inc) 
cU 



< e 



(b) 



< 1. The list of other parameters is given in the caption of 



where (3 = 
point and fi 



2|eSr^|(l 



^(inc)n2wl 



sign(eg™'^'') is the handedness (helicity) of the incident wave. Solutions of 



^cU 



A = arg(teio), Cell = is the eUipticity of the 



equation (68) define the radii of circles containing the C-points. 



Plots depicted in Fig. |4] demonstrate non-monotonic behaviour of the function Re F that 
enter the left-hand side of Eq. (68). 



The experimentally relevant case of the nematic cell filled with the LC mixture E7 is 
shown in Fig.|4]^a). It is characterized by the weak anisotropy of the transmission amplitudes 
\te\ and \to\, so that |1 — |te|/|io|| < 10^"^ for the incidence angles up to 60 degrees. For such 
cells, the phase of the function (69), A, can be approximated by the difference in optical 
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path of the ordinary and extraordinary waves, S = 5e — So- So, the formula 



F(p) ^ Fo(p) = exp(i(5), 5 = Se - 5o, 



(70) 



provides a high accuracy approximation for the expression (69) and, for the LC mixture E7, 
the deviation max \F — Fq\ can be estimated to be better than 5.0 x 10~^. 

The radii of the circles with the C-points. can now be accurately located by solving the 
approximate equation cos 5 = ±P obtained from Eq. (68) by using the expression (70). The 
result can written as two sequences of radii 



\5{po)\ ~ arccos(/?), \^{p^k'^)\ ~ ±arccos(/?) + vr/c, k 



(71) 



where N is determined by the number of solutions and the circles of each sequence, p^^^ and 



Po 



Pl. ', are numbered by the non-negative integer k in non-decreasing order of size, so that 



At p = p\^\ there is a pair of the C-points on each circle of the radius p^' 
ellipticity and the azimuthal angles given by 



< pi-' < pi*' < 



(±) 



with the 



Cell 



(a) 
P=Pk 



-^)^Pinc = Pk, k = 0, 



,iV, 



t* t_ 



(a) ' 



(72) 
(73) 



where pg"'' = po and a G {+, — }. 

The symmetric arrangment of the C-points is a consequence of the symmetry relation 



I^.(P,' 



\^,{p,n + ij)\ 



(74) 



for the amplitudes (58) 



We can now substitute Eq. (58) into the expression for the index (61) and use the rela- 



tions (64) and (68) to recast the result into the form of the index formula 



sign 



dpReFip) 



±(-11 



k 



0,...N, 



(75) 



where the second equality follows because ReF is an oscillating function of p and F{0) = 1. 

Interestingly, Eq. ( 75 ) relates the index of C-points and derivatives of the transmission 
coefficients with respect to the incidence angle {dpf = cos"^ 9incdf /dOinc)- The index is 
determined by the circle number k and alternates in sign starting from Ic = +1/2 which is 
the index of two C-points symetrically arranged in the vicinity of the origin. 

So, the pair of C-points lying on the smallest circle of the radius po can either be of lemon 
or monstar types. Interestingly, when (3 approaches unity and the incident light becomes 
circularly polarized, the radius po vanishes and the C-points merge into one C-point of the 
index +1. 



This is the limiting case where the incident light is circular polarized with e 



(inc) 
ell 



±1 

and a C-point may develop provided that the transmission coefficients meet the isotropy 
condition: tg = ±^o- As it can be seen from Eq. (65) we have = to for the case of normal 
incidence with = sin6'inc = 0. The corresponding C-point is located at the origin and, as 
it was discussed earlier, its index Iq equals +1. 
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At this stage, however, it remains unclear if the above-mentioned C-point is unique. In 
order to clarify this, we shall look more closely at the events that happen when the ellipticity 
parameter (3 defined in Eq. (68) varies from zero to unity. At /3 = 0, the incident light is in 
the state of linear polarization with e^^'^'' = 0, whereas the incident wave is circular polarized 
with le^r^l = 1 at /? = 1. 



The case of linearly polarized incident waves was previously studied in Refs. [291 130] . 
Fig. |6(a)| shows the polarization resolved angular pattern as the polarization ellipse field 
computed at /3 = 0. It can be seen that we have a sequence of + 1 concentric circles each 
containing two pairs of symmetrically arranged C-points. The index and the handedness of 
the C-points alternate in sign along the radial direction. The C-points on the circle of the 
smallest radius, p = po, corresponding to the direction close to the normal incidence are 
found to be of the lemon type with the index equal to +1/2. 

Algebraically, from Eq. (71) it follows that the difference between the radii of the C- 
points p^"*"^ and p^^\, where k = 0,...,A^, vanishes with the ellipticity e^J^^^ approaching 

zero. So, at e^^^^^ = 0, the angular polarization structure is characterized by the presence of 
degeneracy which is a consequence of the additional symmetry relation 



|^,(p,^)| = |*_,(p,7r-^)|. 



(76) 



The symmetry- breaki ng e ffects due to non-vanishing ellipticity of the incident wave are 
illustrated in Figs. 6(b) and 6(c) for small values of This is the case when the rela- 

tion (76) breaks down removing the above degeneracy, so that p^/^^ < p^_,_\ at e^™'^'' ^ 0. 

there is another type of degeneracy which 



As is schematically illustrated by Fig. 5 (a 



takes place when the ellipticity parameter P reaches a local maximum of the function Re F 
such as one marked by box in Fig. |4]^a). In contrast to the above-discussed case of linear 
polarization, the number of C-points changes as the parameter /3 and the ellipticity €^1°^^ 
pass through their critical values, Pc and e^'^-'. 



The picture shown in Fig. 5(a) is typical of a saddle-node (tangent) bifuraction 



l57] . Such bifurcations occur at the incidence angles represented by the values of the radii 
(see Eq. ([56|)) at local extrema of the real part of the function (69). By applying the 



approximation (70), these radii can be estimated from the solutions of the equation 



\6{pi'^)\ ^ TTk, fc = 0. 



N 



(77) 



obtained from Eq. (71) by setting arccos(/5) = 0. 



Given the critical value of the radius, p = p^^\ we can use Eq. (1681) to deduce the relations 



|ReF(pf)) 



(78) 



and to compute the corresponding values of the parameter /3 and the ellipticity: (3 



it' 



and 



^(inc)| 
-cll I 



e'^\ It is also not difficult to show that pg 



and pi") < pt^ 



Our analysis suggests that the ellipse field undergoes a sequence of tangent bifurcations 



as the magnitude of the ellipticity 
of the C-points with p 



pi ^ and p 



,(inc)| 
=ell 



increases from zero to unity. There are two pairs 



p[.^^ that merge and, subsequently, annihilate as the 



governing parameter passes through the bifurcation point with |e 



(inc) I 
ell 



Ac) 



and pl 



(+) 
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In addition, an important consequence of the index formula (75) is that the indicies of 
the C-points involved in a tangent bifurcation differ only in sign, so that the total sum of 
indicies remains intact. Obviously, this result agrees with the conservation law of the total 
topological index. 

Note, however, that, within the approximation (70), the C-points cannot merge and 
annihilate provided the incident wave is not circular polarized. Mathematically, the difficulty 
is that, for the approximate function cos 5, the critical values of the ellipticity are all equal 
to unity, = ei;^ = 1. So, we arrive at the conclusion that the difference 1 — e, is 
determined by the accuracy of the approximation (70). 

For the cell filled with the LC mixture E7, the fisrt three critical values of the ellipticity 
e^r^^ < e''2^ < ef^ turned out to be very close to unity: 1 — ef-* ^ 1.22 x 10~^, 1 — e^/^ ^ 
1.17x 10~^ and 1 — Cg^^ ~ 4.16 x 10~^. The polarization resolved angular patterns presented in 
Fig. 7 are computed at the three different values of e^^^^^: (a) e^^^ < e^^^'^^ = e^/' — 10~^ < e'^^; 

(b) er < 4ir^ = - 10-^ < eS^); (c) e^^^ < e'^.f = e['^ + lO"*' < 1. The pairs of the 
C-points with k = 3, k = 2 and k = 1 are shown to coalesce and disappear in succession as 
the ellipticity passes sequentially through e^f' , e^^ and e^'^''. 

So, in the limit of circular polarization with e^JJ^'^'' = 1, we have the pattern where the 
only C-point is located at the origin. The index of this C-point is Ic = +1. This result 
is in agreement with the Euler-Poincare theorem linking the Euler characteristics, of ^ 
two-dimesional manifold and the total topological charge (index) of a smooth vector field 
on the manifold [8]. In our case, it can be concluded that x = 2/c is the total index, where 
X = +2 is the Euler characteristics of a sphere. 

So far, we considered bifurcations at the local extrema that, according to Eq. (77), are 
in one-to-one correspondence with the extrema of the function ReFo = cos 5. Thus, such 
birfucations are governed by the phase difference between ordinary and extraordinary waves, 
6 = 6e — 6o- As is shown in Fig. [7| they are typically characterized by extremely small 
differences between unity and the critical values of the ellipticity. For this reason they 
cannot be detected experimentally. 

From the above discussion it can be inferred that the critical values might be lowered 
only if the approximation ( [TO]) does not work well. In particular, this is the case for the 
graph of Re F plotted in Fig. |4[b) where the anisotropic material is represented by the KDP 
(potassium dihydrogen phosphate) birefringent crystal. The deviation max \ F — Fq\ can be 

estimated at about 7.0 x 10"^, so that max[l — e^'^^] ~ 0.02. 

The curve shown in Fig. |4]^b) also demonstrates that the breakdown of the approxima- 
tion ( 70 ) may manifest itself in additional extrema formed away from the neighborhood of 
unity. As is illustrated in Fig. 5(b) [ for the region marked by the box in Fig. |4](b), these local 
extrema result in two successive tangent bifurcations leading to creation and annihilation of 
two different pairs of the C-points. 

In contrast to the extrema described by the phase difference S, this is the multireflec- 
tion which is primarily responsible for additional small-scale oscillations. Typically, it is a 
challenging task to resolve accurately these ripple-like and noisy oscillations in polarimetry 
measurements. 
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2. Avoided L-line crossings 



Now we pass on to discussing the L-lines which are the curves of hnear polarization in 



the observation plane. Analytically, substituting Eq. (58) into the relation (63) and using 



the matrix elements (64) give the equation 



2ei;r^ ReF(p) + (1 - [e^Tf ) sin(2^) ImF(p) = 



,(inc)i2^ 
=ell 



(79) 



describing L-lines for the case of the homeotropically aligned NLC cell. 

We begin with the results for the case of linearly polarized incident waves studied in 



Refs. [29l|30]. At e^™'^'* = 0, from Eq. ([79|) there are two straight hues of linear polarization: 



bp'^'^'^ and 



6p"'^^ + 7r/2, where the polarization vectors of incident and transmitted 

waves are parallel, 0^*"^^ = (f)p = (j)p^^\ Other L-lines are circles separating the circles of C- 

points. Interestingly, the radii of the circles, p\^\ can be found as the solutions of Eq. 
with k = l,...N. 

When k is even, it can be concluded that, similar to the straight L-lines, 0p ~ 4>p^^^ 



lik 



is odd, the polarization vector of the transmitted wave rotates with the azimuthal angle of 
the incidence plane and (f)p ~ 0^°'^'' + 20. 
The structure of the L-lines at (3 



is shown in Fig. 6(a) where the L-lines are represented 
by the solid lines bounding the regions of left- and right-handed polarization. As it was 
discussed in the preceding subsection, the case of linear polarization possesses additional 



symmetry described by the relation (76). This symmetry is broken at non-zero values of the 



ellipticity and the case of lin ear polarization is structurally unstable. 



In Fig. 6(b) 



and Fig. 6(c), we show that even small values of the ellipticity, el^^^\ change 



the structure of the intersecting L-lines depicted in Fig. |6 (a)] into a family of non- intersecting 
closed L-lines. The figures clearly demonstrate what might be called, by analogy with the 
well-known avoided level crossings in the quantum mechanics, the effect of avoided L-line 
crossings 

The closed L-lines gradually evolve into the circles as the ellipticity approaches the limit 



(inc) I 
ell 



of circular polarization with |e 

solving the equation: ReF{p) = 0. It follows that p[.~^ at e, 



1. The radii of the C-circles p\^^ can be found by 



,{inc) 
ell 



is equal to p^.^'' at 



1. So, we have a family of the concentric circles describing both the C-points for 



,(inc)| 
Fell I 

linearly polarized incident waves (see Fig. 6(a) ) and the L-lines for circular polarized incident 
waves, (see Fig. [T]). 



IV. RESULTS AND DISCUSSION 

The analytical results presented in the previous section completely characterize the 
polarization resolved angular patterns of light transmitted through the NLC cell for the 
homeotropic orientational structure in terms of the polarization singularities. These pat- 
terns can be conveniently represented as the fields of polarization ellipses and characterize 
the polarization structure behind the conoscopic images measured in experiments with two 
crossed polarizers. 

After computing loci of the C-points and the L-lines, we have deduced a simple formula 



for the index of the C-points ( 75 ) and studied what happen to the polarization singularities 



if the ellipticity and the polarization azimuth of the incident wave vary. It turned out that 
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(a) Experiment at e^"'''' = 0.0 (b) Theory at e^"'"'^ = 0.0 



FIG. 8: (a) Experimentally measured and (b) theoretically computed fields of polarization ellipses 
in the observation plane for the homeotropically aligned cell of the thickness d = 110 /zm filled 
with the NLC mixture E7. The incident wave is linearly polarized with the polarization azimuth 
^(^inc) _ ^ T^liQ C-points of the star, the lemon and the monstar types are marked by stars, 
diamonds and triangles, respectively. L-lines are represented by solid lines. Left-handed and right- 
handed polarization is respectively indicated by solid and open ellipses. 




(a) Experiment at e^^j"''^ = 0.26 (b) Theory at e^j"''^ = 0.26 



FIG. 9: (a) Experimentally measured and (b) theoretically computed fields of polarization ellipses 
for the homeotropically aligned cell filled with the NLC mixture E7. The ellipticity of the incident 
hght is e^j';'''^ = 0.26. 

the C-points may appear and disappear when the pattern undergoes bifurcations at some 
critical values of the ellipticity. In addition, the patterns formed in the neighborhood of zero 
ellipticity (learly polarized incident light) are characterized by the appearance of avoided 
L-line crossings. 
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(a) Experiment at e^)"'^ — 0.65 (b) Theory at e-^' = 0.65 



FIG. 10: (a) Experimentally measured and (b) theoretically computed fields of polarization ellipses 
for the homeotropically aligned cell filled with the NLC mixture E7. The ellipticity of the incident 

(inc) n ctr 

wave IS e„„ = (J. bo. 



Now we compare the experimentally measured and the theoretically calculated polariza- 
tion resolved angular patterns for the homeotropically oriented NLC cell. These patterns 
are shown in Figs. [8 -10 as the fields of polarization ellipses in the observation plane where 
the polar coordinates are defined in Eq. (56) and the scale factor r is taken to be 5 mm. 
The parameters used in our calculations are described in Sec. [TT} They are also listed in the 
caption of Fig. [6] 



In Fig. 



we present the results obtained for the linearly polarized wave with e 



(inc) 
ell 



whereas the cases of elliptically polarized incident hght with e^™^^ = 0.26 and Sl^"^^ 



-- 0.0, 
0.65 



-ell — "-^^ a^^^ 

are shown in Fig. |9] and Fig. [10} respectively. As is evident from the figures, the C-points are 
arranged in chains formed by four rays along which they alternate in sign of the handedness 
and of the index. The L-lines are typically represented by the closed curves separating the 
regions of different polarization handedness. It is also clear that the predictions of the theory 
discussed in Sec. Ill C are in good agreement with the experimental data. 

In addition, by making a comparison between the patterns with different values of el- 
lipticity we arrive at the conclusion that, when the ellipticity changes, the C-points mainly 
move along the radial direction and their azimuthal angles remain approximately constant. 
Theoretically, it can be shown that, at ^p™'^^ = 0, the azimuthal angle l^l^^l given by Eq. (73) 

does not depart appreciably from 7r/4 provided the magnitude of the ellipticity |eg^'^''| is not 
in the immediate vicinity of unity. 

By contrast, as is shown in Fig. [Tj, annihilation of the C-points is accompanied by drastic 
changes in the azimuth of the C-points, whereas the radius is very close to its critical value 
defined in Eq. (77). 

The geometric rearrangements and transformations studied in this paper might be re- 
garded as pseudo-dynamics of the polarization ellipse fields with the ellipticity served as 
a "time" parameter. In three-dimensional space, phase singularity lines, whose geometry 
depends on a "time" parameter were recently investigated in Ref. ^58j 
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In conclusion, we note that the homeotropic ahgnment presents the simplest case of 
anisotropic orientational structure. Our theoretical analysis can be extended to more com- 
plicated cases involving both uniaxial and biaxial anisotropics by using the analytical results 
presented in Appendix |Aj These results will be published elsewhere. 
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APPENDIX A: OPERATOR OF EVOLUTION 



The equations (36) and (38) express the transmission and reflection matrices in terms of 



the linking matrix (34) which is determined by the evolution operator (26). In the case of 



uniform anisotropy, the operator can be computed as the evolution matrix (39) expressed 



in terms of the eigenvalues and the eigenvectors (40) of the matrix (25). 



In this section we consider how the eigenvalue problem can be solved using the basis 
of eigenmodes (normal modes) in a uniformly anisotropic medium. These eigenmodes are 
known to be linearly polarized plane waves characterized by the k dependent refractive 
indices |ll[52l|59]. 



1. Orthogonality relations 



From Eqs. (25a)- (25d), if the dielectric tensor £ is symmetric, eij = e^, we deduce the 



relations for the block matrices 

= Mi2, = M21, = M 

giving the symmetry identity 
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= G M G, G 



I2 
I2 



(Al) 



(A2) 



for the matrix M. 

The left eigenvectors, {V^, Vg, V3, V4}, of the matrix M can be defined as the right 
eigenvectors of its transpose M"^ 



■ = XiY'i ^ ■ V = V ■ A, 



(A3) 



where V' = (V^ V3 V4) is the matrix of the left eigenvectors. It is not difficult to prove 
that the left and right eigenvector form a biorthogonal set 

(V:, V,) = [V,f • V, = 5,, AT, ^ V'^ ■ V = N, (A4) 

where N = diag(A^i, A^2, A^s, ^"4)- 

Eq. (A2) and the definition (A3) can now be combined to yield the relation 

V = G ■ V (A5) 
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linking the matrices of the left and right eigenvectors, V and V. It remains to notice that 

(A6) 



the identities (A4) and (A5) give both the orthogonality relation 

■ G ■ V = N = diag(Ari, N2, N^, N4) 



and the formula for inverse of the eigenvector matrix 



(A7) 



The relations given in Eq. (31) and Eq. (32) represent the special case of isotropic media 
and, thus, immediately follow from Eqs. (A6) and (A7). 

By combining the identities (32) and (A7) with the expression for the evolution opera- 
tor (39), we derive the linking matrix (34) for the case of uniform anisotropy in the following 

W = V^^ ■ exp{-iMh} ■ = N^i . W 



form: 



Wii W12 
-W21 -W22 



w = ■ Wd ■ V, 

V = ■ G ■ V,n, Wd = N-^ ■ exp{-iAh}. 



(A8) 

(A9) 
(AlO) 
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From Eq. ( A8 ) and Eq. (A9 ) the matrices W and Wjj are symmetric, whereas W 
So, we conclude that the transmission matrix T = is symmetric too. 

Note that, for a non-absorbing medium characterized by a real matrix M, ImM 
inverse of the matrix (A8) can be derived by complex conjugation: W^^ = W 
asterisk indicates complex conjugation. In particular, the conservation law 

T+ T + R+ R = I2, 



W 



12- 



0, 

where an 



(All) 



where the superscript f stands for Hermitian conjugation, can be immediately deduced from 
the relation 



W*,-Wn + W*2.W2i 



W+ ■ Wn - W2+ ■ W21 



I9 



(A12) 



2. Uniformly biaxial media 



For plane waves, owing to the Maxwell equation (8b), the electric displacement field is 
transverse, (D, k) = 0. So, assuming that the wave vector k lies in the plane of incidence 
(the x-z plane), the vector D can be conveniently defined by its components in the basis 
{ei(k),e2(k),eo(k) = k} 



B = Diei{k) + D2 e2(k), 

k = eo(k) = k~^k = q'^iq^^ + g^z), 

where ei(k) = g^-'^(g^x — g^jz), e2(k) = y and q = k/k^g^c = n is the refractive index. 



(A13) 
(A14) 



Given the electric displacement field (A13), the electric field can be found from the 
constitutive relation 



E = fir]B, 

77 = ?73 13 + Ar]i d ® d + Ai]2 m ® m. 



(A15) 
(A16) 
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where rji = (/lej) ^, Ar]i = rji — rjs and firj is the inverse dielectric tensor [jit] ■ e = I3). 
The Maxwell equations can now be combined with the relation ( A15 ) to yield the equation 



for the displacement vector D in the form of an eigenvalue problem 

77, -0 = ^-20, r7, = P(k)-r7-P(k), 



(A17) 



where P(k) = I3 — k k is the projector onto the plane normal to the wave vector k. By 



using the expression for the inverse dielectric tensor (A16) we can rewrite the equation (A17) 
in the explicit matrix form: 



[a„ 0-3 + 60-1 - CI2] = 0> 

2azp = Arji{d\ =F d\) + Arj2{fh\ =F m: 

h = Ar]idid2 + Ar]2mim2 
,2 



2J5 



C 



1 

1 



(A18) 

(A19) 
(A20) 
(A21) 



where di = g(d,ei(k)) and nii = g(m, ei(k)) . Then the dispersion relation (Fresnel's equa- 
tion) 



1 - Vsq' 



,2^ 



^V2 [q^ - (rii, q) ) - Arji Ar]2 (d x m, q) =0 



(A22) 



can be derived as the condition for the system of linear equations (A18) to have a non- 
vanishing solution. 

The Fresnel equation describes the wave surface. In our case solving the algebraic equa- 



tion (A22) at q,j^ = n^smOi^c gives the values of the z component of the vector q, q^. 

Generally, there are four roots of Eq. (A22) {qz^\q''^\q^^\qz^^} that form a set of the 



eigenvalues {Ai, A2, A3, A4} defined in Eq. (40). Each root gi*"* corresponds to the eigenwave 
propagating inside the cell with the dimensionless wave vector q^ = (g^, 0, gi*"*) and the 
refractive index Ui = qi. The corresponding polarization vector of the electric displacement 
field is given by 



where 



cos0iei(kj) + sin (pi e2{ki), 
- sin (f)i ei (kj) + cos 0^ 62 (kj 



c (i) > 0, 
c| (i) < 0, 



(A23) 



arg 



= arg(a_ + ib)\ w = 

A?7i((ii + id2y + Ar]2{rhi + im2) 



(i)' 



(A24) 



From Eq. ( A24 ) it is clear that the azimuthal angle 0j becomes indeterminate in the degen- 



erate case when the coefficients a_ and b are both identically equal to zero. Typically, as 
far as the eigenmodes are concerned, this case does not present any fundamental difficulties. 
It just means that the azimuthal angles of the degenerate eigenmodes can be prescribed 
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arbitrarily. Such freedom of choice, however, does not affect the evolution operator which 
remains uniquely defined. 

The procedure to determine the characteristics of the eigenmodes involves the following 
steps: (a) evaluation of the eigenvalues q^f* by solving the Fresnel equation (A22); (b) cal- 
culation of the polarization vectors of the electric displacement field D*-*^ by using the for- 
mula (A23); (c) computing the polarization vectors of the electric and magnetic fields from 
the relations: E^*^ 



/ir7 ■ DW (see Eq. (|Al5|) and /xH^ = x (see Eq. 



As a result, we obtain the eigenvectors expressed as follows 



V,, oc I ^f,) I , 



H 



(i) 



(A25) 



Qz D\ 




Ar7i(q„d)(d,D«) + 



Ar/2(qi,m) (m,D«) 



(A26) 



3. Uniaxial anisotropy 



Now we apply the above procedure to the limiting case of uniaxial anisotropy with t2 = €3 
and Ae2 = A772 = 0. At A?72 = 0, the Fresnel equation (A22) takes the factorized form and 



the values of can be found as roots of two quadratic equations. 

The first equation 1 — rj^q"^ = represents the spherical wave surface. The corresponding 
eigenmodes are known as the ordinary waves. There are two values of 



5; 



(±0) 



± 



nl - ql, 



(A27) 



where n\ = /lea = ^e±, that are equal in value but opposite in sign. When, similar to the 
incident and transmitted waves, the z component of the wave vector (and the vector q) is 
positive, the eigenmode might be called the refracted (forward) eigenwave. In the opposite 
case where, similar to the reflected wave, gi"^ is negative, the eigenmode will be referred to 
as the reflected (backward) eigenwave. So, Eq. (A27) describes two ordinary eigenmodes: 
the refracted eigenwave with q^ = gi^"^ > and the reflected eigenwave with q. 
The second equation 



qi-"^ < 0. 



q^ + Ua{q,dY -n?, = 0, 



(A28) 



where n?, = /lei = fie\\ and Ua = —Arji/rii 



in 



n 



)/n\ is the anisotropy parameter, gives 



the values of for the eigenmodes known as the extraordinary waves. These are given by 



D 



(±e) 



[1 + Uad^^ ^{-Uad^d^q^ ± y/T)] 
'{l + Uad^^ -ql[l + Ua{dl + (i 



(A29) 
(A30) 



where dx = (d, x) and d^ = (d, z) . 

At Ma > {ua < 0), in the x-z plane, Eq. (A28) describes the ellipse with the major 
(minor) semi-axis of the length ny oriented perpendicular to the projection of the director ^ 
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on the plane of incidence {dx,0,dz). The length of minor (major) semi-axis, n± = [nj_ — 
Ua{dy/n\i)'^]~^^'^, depends on the y component of the director and varies from n± to n\\ as dy 
increases from zero to unity. Clearly, degeneracy in refractive indices with Hq = n±e may 
occur only if the director is in the incidence plane (0d = 0). Additionally, the matching 
condition for the x components of q and the director = rim sin 6'inc = ^ngdx = ±^o sin 6^ 
needs to be met. 

The wave vectors and the refractive indices of the normal modes are determined by the 
relation 

q±« = K^c^±a = Qx^ + gi"^"-* z = n±a^±a, a e {o, e}, (A31) 

where n±a = q±a, n±o = Uo = n± is the ordinary refractive index and n±e is the refractive 
index of the extraordinary wave propagating along the unit vector k-te. 

By substituting Ari2 = into Eq. (A23) we can obtain the polarization vectors of the 
electric displacement field |29l |30] : 



D(±°) = -ql,q±o X d = -nlq±o X d, (A32a) 

D(±^) = gLP(k±e) ■ d = g|,d - (d,q±e)q±e. (A32b) 

Then, following the procedure described at the end of the preceding section, we find the 
polarization vectors of the electric field for the eigenmodes 

^-1e(±-) = d -nl2(d,q±,)q±,. (A33b) 
The result for the magnetic field of the normal modes is 

H(±°) = < d - (d, q±,) q±o, (A34a) 

H(±^) = q±e X d. (A34b) 

Solution of the eigenvalue problem can now be presented in the matrix form: 

A = diag(g(+'=), g(-), V = (^^+ , (A35) 



where 



E± = I (±e) J±o) h H± = I . (A36) 



Owing to the orthogonality relations (A6) the matrix 

N = ■ G ■ V = diag(N+, N_) (A37) 



is diagonal and its non-vanishing elements 

N± = {^^^ J , = 2(z,E(") X hW) = 2(Ei,"\H?)) 
are proportional to the normal components of the Poynting vector of the eigenmodes. 



(A38) 
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Substituting the expressions for the electric and magnetic fields of the eigenmodes [see 
Eq. (A33) and Eq. (A34), respectively] into Eq. (ASS) gives the diagonal elements of the 
matrices N. and N_ 



/.-^ N^o = Iqf^^ {n\ - (q±o, d) , (A39) 
/i-^ iV±, = 2n;2{#^) {n\ - (q±„ d)'] + 4(q±e, d) [n^ - n\\\ (A40) 

The eigenvalues of ordinary and extraordinary waves that enter the matrix of eigenvalues, 
A, are given by Eq. (A27) and Eq. (A29), respectively. The elements of the 2x2 block 
matrices (A36) that define the matrix of eigenvectors (A35) are the lateral components of 
the vectors (|A33l) and (|A34l). 



These analytical expressions can now be substituted into the linking matrix (AS) to write 
explicitly the elements of the transfer matrix (36 ). For the homeotropic structure with d = z, 
it is rather straightforward to show that the result is given by Eqs. (41)-(46). 
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